Is diffusion anomalous in two-dimensional Yukawa liquids? 
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There have recently been many predictions of "superdiffusion" in two-dimensional strongly cou- 
pled Yukawa systems, both by computer simulations and in dusty plasma experiments, with substan- 
tially varying diffusion exponents. Here we show that the results crucially depend on the strength of 
dissipation and the time instant of the measurement. For sufhciently large friction even subdiffusion 
is possible. Ifowever, there are strong indications that, in the long-time limit, anomalous diffusion 
vanishes for dissipative as well as for frictionless systems. 
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Two-dimensional (2D) systems are important models 
for a large variety of processes on liquid or solid surfaces 
as well as for low-dimensional structures in condensed 
matter. It is long known that reduction of the dimension- 
ality to two is accompanied by a number of anomalies in 
the collective properties of many-body systems including 
Bose condensation and superfluidity, the quantum Hall 
effect or phase transitions (Kosterlitz-Thouless scenario) . 
Even purely classical systems have been predicted to ex- 
hibit anomalies in 2D which are seen e.g. in fundamental 
transport properties such as diffusion ^1 . 

Among classical 2D systems, monolayers of particles 
interacting via a Coulomb or Yukawa potential which 
can be produced in ion traps, e.g. ^ and dusty plasmas 
[3] constitute a particularly interesting generic example 
which allows to study 2D anomalies with unprecedented 
accuracy. It is, therefore, not surprising that diffusion 
processes in 2D Yukawa liquids (2DYL) have, over the 
last decade, attracted considerable interest, both in the- 
oretical works [3 El |71 IHJ 13 Eni HI] and experiments 
in dusty plasmas |4l|7l[l2l[13l[l4l[l5l[16l[17l[18l[19]. 
While some authors found no peculiarities [U [S], the 
majority of these works has reported significant devia- 
tions from "normal" diffusion, i.e. from Pick's law and 
Einstein's formula for the mean-squared displacement 
(MSD), Ur ~ with a = 1. There were many obser- 
vations of superdiffusion {a > 1), i.e. enhanced diffu- 
sion [iHiHiiiiiniinjiiaiiiiisiiiniiniiTsiiTS], but 

in some cases also of reduced diffusion, i.e. subdiffusion 
{a < 1) [51 [5]. The strength of the diffusion can be 
conveniently quantified by the value of the exponent a. 
However, there is a substantial scatter of the reported ex- 
perimental and theoretical data, ranging from below 1.0 
to values as high as 1.3. The origin of these differences is 
unknown, and even the existence of anomalous diffusion 
is being debated. 

The theoretical investigations were based on equilib- 
rium computer simulations and have concentrated on the 
idealized 2D as well as on the quasi-2D case [8|, both with 
and without dissipation included. Possible sources of de- 
viations lie in differences in the complex plasma condi- 



tions such as screening strength k, coupling parameter P 
and neutral gas friction. In dusty plasma experiments [S], 
moreover, various forms of friction are always present as 
are additional energy sources arising from the plasma dis- 
charge giving rise to a nonequilibrium driven-dissipative 
system. 

The aim of this Letter is to resolve the puzzle about the 
character of diffusion in 2DYL. To this end, we systemat- 
ically analyze the effects of the interaction range (screen- 
ing), coupling (temperature) and dissipation. Purther- 
more, it is well known that the diffusion exponent a varies 
in time: initially it equals 2 (free ballistic motion) and 
lateron it is expected to approach some fixed value which 
determines whether the system exhibits normal, sub- or 
superdiffusion. It has been noted UHl HH HI HI US] 
that this value depends on the choice of the time win- 
dow during which it is recorded. This may be another 
explanation for the large scatter of a values reported pre- 
viously but raises the question about the correct proce- 
dure. We therefore, explore the whole time dependence of 
a in detail extending the simulations to very long times. 
The results can be summarized as follows: the charac- 
ter of diffusion depends on the magnitude of dissipation. 
While for weak friction, superdiffusion is observed, an in- 
crease of friction eventually gives rise to normal diffusion 
and subdiffusion. However, this turns out to be only a 
transient phenomenon. Por sufficiently large observation 
times, diffusion transforms to normal diffusion. 

Model and simulation idea. Our 2DYL consists of 
N particles contained in a quadratic monolayer with pe- 
riodic boundary conditions. The particle motion is mod- 
elled by coupled Langevin equations, 

niifi ^ Fi- m^Dvi + i^ I. ..N (1) 

where 9 is the friction coefficient and Fi the Yukawa force, 

47r£o^V r J 

with the Debye screening length Xd and charge Q. yi{t) is 
a Gaussian white noise with zero mean and the standard 
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FIG. 1; (Color online) Diffusion exponent a as a function of 
r'^'^' for different values of k (1.0,2.0,2.5,3) obtained from a 
best fit of Eq. ^ in the time window ujpt G [100, 320] The 
three ticks on the upper x-axis indicate the three values of 
r"'=' used in Figs. 

deviation {ya,i{to)ypj{to+t)) = 2kBTV5ij5af3 5{t), where 
a, /? € {x, y} and T is the temperature. We choose as the 
unit of length the Wigner-Seitz radius a^s = (n7r)^^/^, 
where n is the areal density, and as the unit of time the 
inverse of the plasma frequency ujp = {Q"^ /iireoma^gY/'^ . 
In thermodynamic equilibrium, the system is thus fully 
described by three parameters ~ the dimensionless in- 
verse screening length k — Qws/^d, the friction coeffi- 
cient v = P/ujp and the Coulomb coupling parameter 
r = Q^/{4:TreoakBT). Eq. Q is solved by standard 
Langevin dynamics |20j up to a maximum observation 
time tabs which is limited by the condition that no col- 
lective oscillations, e.g., sound waves, should be able to 
traverse the entire simulation box of length L during the 
measurement [11]. Thus, tots < L/vg should be satisfied 
where Vg is the (sound) velocity of the fastest mode which 
for 2DYL is well known [1TJ[21|13]. By simulating large 
systems with 100,000 particles, we achieve tobs = 2,500 
plasma cycles, extending previous results [51 HI IHl HH] by 
one order of magnitude. 

Mean squared displacement. In examining the 
character of diffusion, the most important quantity is the 
MSD, Ur{t) — {\f{t) — r(to)P)Ar, where the averaging is 
over all particles. A general parametrization is 

Ur{t) ^ Doe^''> , (2) 

where in case of normal diffusion, at long times a ap- 
proaches unity, giving rise to the diffusion coefficient 
D = Ur{t)/4:t. In contrast, in the case of anomalous dif- 
fusion, a differs from unity: a > 1 (a < 1) is associated 
with superdiffusion (subdiffusion) . 

In a first series of simulations with fewer particles and 
zero friction, we establish the overall diffusion trends in 
2DYL, i.e., the dependence of a on both the coupling 
strength and the range 1 /k of the pair interaction. To be 
able to compare the data for different k, we introduce the 
relative coupling F"'"' by normalizing the Coulomb cou- 



pling parameter to the crystallization point Tc{k), i.e., 
prci _ Y/Tc [10]. This is only one possible definition of 
the physical coupling - see e.g. Refs. [22l[24l[25] for other 
definitions - but a particularly physically intuitive one. 
The results of these simulations are shown in Fig. [T] A 
first observation is that a is larger than unity for a large 
range of F'''^\ clearly indicating superdiffusion confirming 
earlier simulations [6] [9l [10]. Our data are significantly 
more comprehensive covering the full range of coupling 
strengths from the gas phase to the strongly coupled 
liquid. The dependence is non-monotonic with a max- 
imum around F'''^' — 0.2. Our data show that superdiffu- 
sion is strongly dependent on the degree of correlations 
in the system. For weak coupling, the system is domi- 
nated by binary interactions, and a is only slightly larger 
than unity. When the coupling is increased, superdiffu- 
sion becomes more pronounced because collective effects 
grow. After a maximum, the value of a is again decreas- 
ing since particle movement is increasingly hindered by 
entrapment in local potential minima ( "cages" ) and the 
onset of crystallization. A second conclusion from Fig. [T] 
is the weak but systematic dependence of a on k in the 
range k = 1 ... 3. This trend has been first observed in 
[To] for three values of F''''' and is here confirmed for the 
whole range of coupling strengths. 

Effect of dissipation. We now turn to the infiuence 
of friction and to the long-time behavior of diffusion. To 
achieve long simulation times, tobs, it is advantageous 
to consider large n values since the sound speed Vg di- 
minishes rapidly with k. Below we will, therefore, con- 
centrate on the case k — 3 which - due to the weak k 
dependence observed in Fig. [T] - is representative for a 
2DYL. This allows us to perform extensive large simula- 
tions with TV 100, 000 up to twp = 2, 500. 

In Fig. [2|a), we show Ur{t) for three values of the fric- 
tion coefficient z/ at a fixed coupling F''°' = 0.17 close to 
the maximum of the curves in Fig. [l] corresponding to 
a moderately coupled liquid state. The ballistic regime, 
LUpt ^ 5, is followed by a transition period after which 
the MSD eventually seems to approach an asymptotic 
behavior characterized by different slopes for different v. 
From these curves it is clear that the common method to 
extract a from the average slope of the MSD in a given 
time window, e.g., [5] [HI US], can be rather ambiguous, 
depending on the width and position of the latter. There- 
fore, here we avoid any averaging and, instead, study an 
instantaneous diffusion exponent a'"*'' as the derivative 
a'"«'(i) = 9i„tlnK(lnt)]. 

Fig. ^h) shows a""*(t) for the same curves (we drop 
the superscript in the following). Evidently, the instanta- 
neous exponent a(t) is much more sensitive to the time- 
dependence of the diffusive motion which turns out to be 
non-trivial. Consider first the curve for the friction value 
I' = 0.08. Here, a{t), after a transient period of about 20 
plasma periods, approaches its asymptotic value a = 1, 
i.e. the system exhibits normal diffusion. In contrast, for 
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FIG, 2: (Color online) (a) Urit) for T''^' 
0,001,0,08,1,0, (b) Corresponding 



0,17 and u 



the low-friction case, v = 0,001, after a transient of sim- 
ilar length a{t) remains significantly larger than unity 
corresponding to superdifFusive behavior. However, no 
constant value a(t) is reached within our simulation time, 
rather a{t) continues to decrease monotonically towards 
unity. Even more interesting is the behavior at large 
friction, ly — 1: after the ballistic regime, a broad in- 
termediate phase is observed where the motion changes 
from superdifFusion to subdifFusion, The latter persist 
for 10 < ujpt < 2000 after which the diffusion is normal 
(Here, and in the following, we use the threshold a = 1,03 
[a = 0,97] to locate the boundary between normal diffu- 
sion and superdifFusion [subdifFusion]), 

Transient and long-time behavior for different 
couplings. To systematically examine the influence of 
dissipation on a(t), we repeat these simulations for a 
larger number of friction coefficients and different cou- 
plings. The entire time dependence of a(t) for arbitrary 
v at r''"^ = 0.17 is comprised in Fig,|3] A first observation 
is that, in all cases, dissipation ultimately induces a tran- 
sition to normal diffusion, Superdiffusion may exist, but 
it is a transient effect which is possible only at sufficiently 
low friction. The time window of superdifFusion rapidly 
decreases as dissipation grows. At high friction, instead, 
transient subdifFusion is observed within a time window 
which is growing with increasing friction. Only within a 
relatively small friction interval a direct crossover from 
ballistic motion to normal diffusion is observed. This is 
an interesting special case where the combined effect of 
friction and thermal fluctuations prevents the build-up 
of collective motions in the liquid. The corresponding 
graphs for the limits of very strong and very weak cou- 
pling, V"^^ = 0,43 and 0,043, respectively, are shown in 
Figs, [4] Here, the overall picture is the same as for 
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FIG, 3: (Color online) Instantaneous diffusion exponent 
a{t) as a function of friction coefficient 1/ and time cupt for 
F'^''^ = 0,17, (a) Density representation with contour lines for 
a = 1,50, 1,03, 0,97 (b) Schematic representation with labelled 
regions of super-, sub- and normal diffusion. 



prci _ Q j^y with the main modification that, for lower 
[stronger] coupling, the region of transient superdifFusiv- 
ity is extended towards higher [lower] friction, while the 
subdifFusive region is diminished [increased]. Thus we 
conclude that we have established the general trends of 
time-dependent diffusion processes in dissipative 2DYL, 
Now the interesting question arises whether the emer- 
gence of normal diffusion at very long time scales is the 
sole consequence of dissipation or if this is an intrinsic 
property of 2DYL which is observed also in the friction- 
less limit. To this end, we have performed an additional 
series of simulations for v = and a broad range of cou- 
plings the results of which are displayed in Fig, |5] For 
high couplings, a{t) shows oscillations during the tran- 
sient regime with a frequency close to the Einstein fre- 
quency [24] which are due to caged motion of the parti- 
cles. This transient regime is followed by a long period 
of superdiffusive motion which seems to tend towards 
normal diffusion for longer times. On the other hand, 
at weak coupling, a(t) does not exhibit oscillations but 
rather decays monotonically from 2 to 1, Whether in 
the long time limit normal diffusion will emerge, can- 
not be answered definitely from our data due to the lim- 
ited observation time. However, an extrapolation of the 
data seems to indicate that a(t) converges towards unity 
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FIG. 4: (Color online) Instantaneous diffusion exponent a{t) 
as a function of friction coefficient ly and time ujpt for F"^"^' = 
0.43 (top) and F''^' = 0.043 (bottom). 



for t > 10^0-"^ 



for 



within two more time decades, i.e. 
all coupling strengths. 

In summary, in this Letter we have presented a com- 
prehensive numerical analysis of diffusion in 2DYL, both 
with and without friction. We have clarified the effects 
of coupling, interaction range and dissipation and, most 
importantly, the time dependence of the diffusion expo- 
nent over long time scales. Our results explain the partly 
contradicting results of earlier studies and confirm the ob- 
servation of superdiffusive behavior. However, this turns 
out to be a transient effect existing only in a finite time 
window and below a critical dissipation which depends 
on the coupling strength. For larger dissipation, instead, 
transient subdiffusion is observed. For the limit of non- 
dissipative 2DYL our simulation results provide strong 
indications that here anomalous diffusion vanishes as well 
on sufficiently long time scales ujpt > 10^. We underline 
that our results do not rule out superdiffusion in dusty 
plasma experiments because they are performed under 
stationary nonequilibrium conditions. To verify the rel- 
evance of nonequilibrium states, we suggest to perform 
measurements of the time dependence of a and compare 
it to our data. 

The observed nontrivial anomalous diffusion of clas- 



sical particles arises from their confinement to two di- 
mensions. Similar transient anomalies should show up 
in other transport properties. Moreover, it may be ex- 
pected that 2D anomalies of correlated quantum systems 
will exhibit similar time-dependent modifications as well 
which should influence their excitation spectra. 




FIG. 5: (Color online) Instantaneous diffusion exponent a{t) 
for u = 0. Subsequent curves differ in F''^' by a factor of 1.41. 
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